Side-jumps in the spin-Hall effect: construction of the Boltzmann collision integral 



o 
o 

(N 
> 

o 
O 
in 



Dimitrie Culcer,^ E. M. Hankiewicz,^ Giovanni Vignale,'^ and R. Winkler''^ ^ 

'Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park MD20742-4111 
^Institut fiir Theoretische Physik und Astrophysik, 
Universitdt Wiirzburg, 97074 Wiirzburg, Germany 
^Department of Physics and Astronomy, University of Missouri, Columbia, Missouri 65211 
^Materials Science Division, Argonne National Laboratory, Argonne, IL 60439 
''Northern Illinois University, De Kalb, Illinois 60115 
(Dated: October 15, 2009) 

We present a systematic derivation of the side-jump contribution to the spin-Hall current in 
systems without band structure spin-orbit interactions, focusing on the construction of the coUision 
integral for the Boltzmann equation. Starting from the quantum Liouville equation for the density 
operator we derive an equation describing the dynamics of the density matrix in the first Born 
approximation and to first order in the driving electric field. Elastic scattering requires conservation 
of the total energy, including the spin-orbit interaction energy with the electric field: this results 
in a first correction to the customary collision integral found in the Born approximation. A second 
correction is due to the change in the carrier position during coUisions. It stems from the part of the 
density matrix off-diagonal in wave vector. The two corrections to the collision integral add up and 
are responsible for the total side-jump contribution to the spin-Hall current. The spin-orbit-induced 
correction to the velocity operator also contains terms diagonal and off-diagonal in momentum space, 
which together involve the total force acting on the system. This force is explicitly shown to vanish 
(on the average) in the steady state: thus the total contribution to the spin-Hall current due to the 
additional terms in the velocity operator is zero. 
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I. INTRODUCTION 

Semiconductor spin electronics is an active area of re- 
search in which both theory and experiment have made 
substantial progress during the past decade. ^'^ The re- 
cent focus on electrically-induced phenomena in sys- 
tems with spin-orbit interactions has brought to light 
many unexplored and fascinating facets of semiconduc- 
tor physics. Considerable progress has been made in past 
years in the electrical manipulation of spins in semicon- 
ductors, where experimental and theoretical work have 
yielded the prediction"^"^ and discovery of the spin-Hall 
effect. ^""^^ The spin-Hall effect, which is the main fo- 
cus of this paper, is the generation of a transverse spin- 
current^'^"* at the edges of the sample as a response to 
an external electric field. Such a spin current leads to a 
spin-accumulation at the edges of the sample, and the re- 
lationship between spin currents and spin accumulation 
is nontrivial.*^'*^ The first observations of the spin-Hall 
effect were followed by the measurement of the spin-Hall 
effect at room temperature by optical techniques*^ and 
the first successful measurement of the spin-Hall effect in 
transport in balhstic HgTe/HgCdTe quantum wells.** 

Two main mechanisms have been shown to be re- 
sponsible for the spin-Hall effect. The presence of spin- 
orbit coupling in the impurity potentials gives rise to 
contributions to the spin-Hall effect which are termed 
extrinsic.^^^'^^'^'^^^^ The spin-Hall effect observed in Refs. 
10,11 was shown by Engel et al.,^^ Tse and Das Sarma^° 
and Hankiewicz and Vignale^^ to be due to extrinsic 
mechanisms. Band structure spin-orbit interactions yield 
contributions termed intrinsic.^''' '"^^^"^^ The spin-Hall ef- 



fect observed in Ref. 12,18 is believed to be due to in- 
trinsic mechanisms. Extrinsic and intrinsic mechanisms 
were broadly discussed in the context of the anomalous 
Hall effect (AHE), which is the generation of a transverse 
charge and spin-polarization current in response to an 
electric field in a ferromagnetic medium. ^^"'*'*' In fact, ex- 
trinsic contributions to the anomalous Hall and spin-Hall 
effects are closely related. ^^'^^ The interplay of intrinsic 
and extrinsic contributions is a complicated problem. It 
was first addressed by Tse and Das Sarma'*^ using an ap- 
proach based on the diagrammatic Kubo formula. This 
was followed by a series of publications'*^"^* and this topic 
continues to be an active area of research. However, our 
focus in this work is on the case of extrinsic spin-orbit 
interactions alone and specifically the way they may be 
obtained from a kinetic equation approach. 

Extrinsic contributions to the spin current are of two 
kinds. The first and more intuitive contribution arises 
from the asymmetric scattering of up and down spins 
known as skew scattering. This effect is found be- 
yond the first Born approximation, i.e., provided one 
goes to at least third order in the electron-impurity 
potential. *^'^°''^^ The effect appears naturally in the 
standard Boltzmann collision integral^^''^^ provided one 
goes beyond the first-Born approximation. The associ- 
ated spin-Hall conductivity scales with disorder as the 
ordinary Drude conductivity*^'^° (i.e., proportional to 
the electron-impurity scattering time), although it is of 
course much smaller. 

The second extrinsic contribution has been known in 
the literature as side-jump and has two main character- 
istics: (i) it appears already in the first Born approxi- 
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mation for the electron-impurity potential, and (ii) the 
associated spin Hall conductivity is independent of the 
electron-impurity scattering time a surprising univer- 
sality which will be fully explained below. In contrast 
to skew-scattering, the side-jump mechanism cannot be 
derived straightforwardly from the standard form of the 
Boltzmann equation. Very early on, Luttinger^° calcu- 
lated the side jump contribution to the charge conductiv- 
ity using a recursive density matrix approach, providing 
a thorough calculation yet leaving many questions unan- 
swered concerning the physical picture behind different 
contributions. Later, Berger''^"^''' used a wave-packet for- 
malism to build a semiclassical picture of side jump, 
identifying it with a shift of the center of mass of the 
wave packet during collisions. Nozieres and Lewiner^^ 
used this picture and carefully studied the side jump 
contributions within a Boltzmann approach. These au- 
thors accounted for six contributions, some of which can- 
cel each other; yet they did not associate any physi- 
cal interpretation with the cancelations. In Ref. 15 the 
skew-scattering and side-jump terms were identified also 
within the framework of a drift-diffusion approach. In 
Rcf. 23 it was claimed that the proper way to describe 
the side-jump effect is to replace the band energy in the 
usual Boltzmann collision integral by a modified band 
energy which includes corrections due to the spin-orbit 
interaction with the electric field. However, the validity 
of this claim has not been formally shown to date. Skew 
scattering and side jump were rigorously studied by Tse 
and Das Sarma^'^ using a diagrammatic Kubo formalism, 
demonstrating that the side-jump term can be identified 
with an anomalous current which gives rise to a renor- 
malization of the current vertex. A derivation based on 
the Kubo formula was also presented in Ref. 49, finding 
that the side-jump contribution for the conduction band 
is independent of disorder and of the Coulomb potential 
to all orders in the strength of these interactions. 

The purpose of this paper is to show that a rigorous 
alternative derivation of the side jump within the frame- 
work of the kinetic equation which does not require intu- 
itive approaches, and follows cleanly and clearly from the 
Boltzmann equation provided one constructs the collision 
integral with care. This construction, however, requires 
that we go beyond the Boltzmann equation formalism 
and resort to the full quantum Liouville equation for the 
one-particle density matrix, which we treat to second or- 
der in the electron-impurity potential and to first order 
in the electric field. In order to focus on the essential 
physics, we assume that band structure spin-orbit inter- 
actions are negligible. Thus we only take into account 
spin-orbit interactions with the impurities and with the 
external electric field that drives the current. 

Starting from the quantum Liouville equation we show 
rigorously that the conservation of energy in the Boltz- 
mann collision integral must be modified to include the 
effective band energy, i.e., the bare band energy plus the 
spin-orbit interaction with the electric field. This modi- 
fication is reflected in the appearance of a correction to 



the scattering term usiially found in the Born approxima- 
tion. This additional term, which is spin-dependent and 
linear in the electric field, acts as a source for the spin 
current, yielding one half of the side jump contribution. 

The other half emerges when one takes into account 
the change in the position of the particle during colli- 
sions with impurities. The semiclassical picture of the 
renormalization of the trajectory of the electron during 
collisions will be derived rigorously from the scattering 
of electrons off the impurity potential which involves the 
terms off-diagonal in the wave vector in the density ma- 
trix. It turns out that this process and the modified 
conservation of energy give equal contributions to the 
spin-Hall current, so that the full side jump contribution 
can indeed be obtained by including twice the spin-orbit 
interaction energy with the electric field in the (5-function 
of conservation of energy in the ordinary Boltzmann col- 
lision integral, as suggested previously in Ref. 23. How- 
ever, the present derivation is rigorous, and constitutes 
formal validation of this heuristic approach. We note 
that the factor of two emerges naturally from the dia- 
grammatic Kubo formula, as demonstrated in Ref. 20, 
when one takes into account the vertex renormalization 
of the spin and charge currents. 

In addition to constructing the picture of side jump 
outlined above, the analysis expounded in this work sup- 
ports the argument that the full velocity operator con- 
tains the net force acting on the system, which must van- 
ish on physical grounds. To provide a formal derivation 
of this fact, we show that the total contribution of the 
additional terms in the velocity operator to the spin-Hall 
current is indeed zero. This is due to the fact that the 
velocity operator contains a spin-dependent term linear 
in the electric field, which reflects the spin-orbit interac- 
tion with the electric held, as well as an additional spin- 
dependent term off-diagonal in wave vector, which re- 
flects the spin-orbit interaction with the impurities. The 
importance of this latter term was recognized already in 
diagrammatic linear response theory in Ref. 20. Within 
the kinetic equation formalism, the two terms in the ve- 
locity operator produce equal and opposite contributions 
to the spin-Hall current which cancel each other out. The 
physical interpretation of this fact is that, in the steady 
state, the average force acting on an electron must van- 
ish. 

In this paper, therefore, we strive to provide an imder- 
standing of the side-jump mechanism in the absence of 
intrinsic spin precession due to band-structure spin-orbit 
coupling, which is rigorous and at the same time physi- 
cal. We believe that a rigorous physical understanding is 
a first step towards building a consistent picture, within 
the kinetic equation framework, of the interplay of spin 
precession due to band structure spin-orbit coupling and 
spin-orbit coupling due to impurities. This interplay has 
been studied by Tse and Das Sarma"'^ based on a dia- 
grammatic linear-response approach, and by Hu et aZ^^ 
numerically. More recently Hankiewicz and Vignale^^ 
constructed a phase diagram of this problem, while Rai- 
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mondi and Schwab'*** employed a Keldysh Green's func- 
tion technique. The long-term aim of this work is to build 
a rigorous understanding, based on the kinetic equation, 
of intrinsic spin precession, skew scattering and side jump 
on an equal footing. 

This article is organized as follows. In section II a ki- 
netic equation is derived starting from the quantum Li- 
ouville equation for the density operator. In section III 
we explicitly construct the collision integral, including 
all contributions in the first Born approximation arising 
from the modification of the position operator. Section 
IV discusses the contributions of the side-jump mecha- 
nism to the spin-Hall effect. It demonstrates that the 
corrections to the velocity operator do not contribute to 
the spin current. We end with conclusions. 



II. TIME EVOLUTION OF THE DENSITY 
MATRIX 

We outline in this section the formalism that will be 
used to determine the collision integral and the way it 
gives rise to the side jump spin-Hall current. The system 
is described by a density operator p, which obeys the 
quantum Liouville equation 



dp 
'di 

In this equation 



-£ + ^[Ho + He + Hu,p]=0. 



21.2 



Ho 



2m* 



(1) 



(2) 



is the Hamiltonian for a parabolic conduction band, with 
m* the carrier effective mass, while 



He = eE -f + XeE ■ & X k 
= eE-r + ^&-Ak 



(3a) 
(3b) 



represents the interaction with a constant uniform exter- 
nal electric field E. The term He includes both the direct 
interaction and the interaction via the (material depen- 
dent) spin-orbit coupling of strength A. We have used 
the notation = 2Xek x E for the effective field char- 
acterizing the spin-dependent part of this interaction, a 
term which, as we shall see later, gives rise to one half of 
the side jump contribution to the spin-Hall conductivity. 
Finally 



Hu = U{r) + XVUir) ■ & x k 



(4) 



denotes the interaction with a set of randomly distributed 
impurities, again both directly and via spin-orbit cou- 
pling. 

We project the Liouville equation onto a set of time- 
independent states {|fes)} of definite wave vector k and 

spin orientation s ~ ±1 along the z-axis. The matrix 
elements of p in this basis form the spin density matrix 



and arc written as pfefc/ = pl^, = (fc,s|/5|fc's'), with cor- 
responding notations for the matrix elements of i/o) He 
and Hij. We assume that impurities are uncorrelated and 
the normalization is such that the configurational average 
of {ks\Hij\k' s')l^k' s'\Hu\ks) is (rii/F) \lJkw^8ss', where 
n, is the impurity density, V the crystal volume, and 



Ukk' = Ukk' (l - «Acr ■ k X k') 



(5) 



Here Ukk' are the matrix elements of the electron- 
impurity potential U{f) between plane waves, while U^k' 
is reserved for the total potential of a single impurity in- 
cluding the spin-orbit contribution. In what follows spin 
indices will be suppressed and all quantities are assumed 
to be matrices in spin space. Note that we use the conven- 
tion that Ukk' and Ukk' have units of energy x volume. 

The density matrix p is divided into a part diago- 
nal in k and a part off-diagonal in fe, given by pkk' = 
fk Skk' + 9kk' ■ These two parts of p satisfy a set of cou- 
pled equations 



[Hu,g] 



kk 



(6a) 



dg, 



kk' 



dt 



kk' 



h 

[Hu,f]kk' - [Hu,9\kk',{Q'^) 



where H = Hq + He- In the first Born approximation 
the solution to Eq. (6b) for gkk' is 

/•OO 

gkk' = -l lini / dt' e-"*' e-^"''/^[Hu, f{t-t')] e^"''/%^, 

(7) 

where ?7 > is a regularization factor. We are considering 
variations which are slow on the scale of the momentum 
relaxation time, thus we approximate f{t — t') in the inte- 
gral by f{t),^° which is written simply as /, and satisfies 
the equation 



Afk) 

where the scattering term J(/fe) is 



J{fk) = {i/K)[Hu,9\kk 



(8) 



(9a) 



1 /■oo 

^ lim dt' e-"*' [Hu, e-'"''l^[Hu, f] e'"'''%^^. 

(9b) 



Equations (8) and (9) describe the dynamics of the den- 
sity matrix and constitute the complete set of tools we 
require in order to derive the kinetic equation satisfied 
by the density matrix in an electric field, including the 
side-jump contribution to the scattering term due to the 
modification of the position operator by the spin-orbit 
interaction. 
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III. DERIVATION OF THE COLLISION 
INTEGRAL 

We want to evaluate further the collision integral (9). 
For this purpose, we decompose the matrix fk into a 
part scalar in spin space and a spin-dependent part, thus 
fk = nfc 1 + Sfe, with 1 the identity matrix and Sk ex- 
pressible in terms of the Pauli matrices. We will show in 
the following that, to first order in A and in the electric 
field, we can write the scattering term as 

J{fk) = ioK) + J^iiuk) + J^j(nfe) + Jo{Sk). (10) 

The first of these terms comes from the band Hamilto- 
nian Hq, is a scalar in spin space and does not depend 
on A or on the electric field. The second term reflects 
the fact that the total energy must be conserved during 
scattering events, including the second term in Eq. (3). 
The third term comes from the direct interaction with 
the electric field eE ■ r, and arises because r fails to com- 
mute with the spin-orbit interaction with the impurities. 
Physically this reflects the change in r during a collision 
with an impurity. The resulting variation of eE ■ r also 
contributes to the overall energy balance. Both J^(nfc) 

and Jgj(nfe) are spin-dependent and are linear in A and 



the electric field E. We will discuss each of the contri- 
butions in turn below. The last scattering term, Jo(S'fc), 
will be important in the kinetic equation below in de- 
termining the steady-state correction linear in E to the 
spin-dependent part of the density matrix, and thus to 
the spin-Hall current. 

A. Scattering correction due to the conservation of 
the modified carrier energy 

In this subsection we focus on the part of the scattering 
term which is linear in the electric field and arises from 
the addition of the spin-orbit interaction with the electric 
field i (T • Afe to the particle energy eofe = ft^fc^/(2m*). 
Since we are working to first order in A and in the elec- 
tric field, in this subsection we only need to consider the 
scalar part of the impurity potential, Ukk'- Moreover 
the time evolution operators in this subsection include 
the side jump energy ^ cr • A^, but not the term eE ■ r, 
which will be considered in the next subsection. The ma- 
trix elements of H in the exponents of the time evolution 
operators are thus diagonal in k. 

Writing out all the terms in the double commutator 
(9b) we find 



1 [Hu, e-'"'''\Hu, f] e'^''/%k = ^ (Ukk'S-^'''^"' ^^Uk'kfke'"'''' - Ukk'e-^"^"' ^''fk'Uk'k e^^'=*'/'^ 

-e-'"-''/''Ukk'fk'e'"'''*'/''llk'k + e-'"''''"'fkUkk'e'"'''*'"'Uk^k) , 

(11) 

where Hk = eok + ^ • and d is the dimensionality of the system. By expanding the exponentials of Pauli matrices, 
the product of time evolution operators can be written, to flrst order in A, as 



g-iHk't' /h^iHkt' /h _ gi(eofc-eofc')*7?i 



Akt' A'j^t' . ^ Akt' . A'^t' , . ^ . Akt' A'^t' 
cos — — — COS tcr ■ Afc/ cos — — — sin —ft — hicr ■ At sm — — — cos 



2h 2h 2h 2h " 2h 2h 

(12) 

where A^ is a unit vector in Afc direction. 

The only task that remains is integration over the time variable t', giving a series of (5- functions reflecting energy 
conservation. The overall result for this scattering term, to first order in A, can be decomposed into a scalar part 
•^0 {nk) independent of A, and spin-dependent parts J^- {uk) + Jq (Sk)- These parts may be written as follows 

Jo (rife) = ^ J-0^\Ukk'\\nk-nk')[5{e+-e'+)+6{e.-e'_) + 6{e+-e'_) + 5{e.-e'+)] (13a) 
J^iink) = ^ |^|WfefeH^(nfe-nfeO{eT.(Afe + AfeO[<5(e+-eV)-^(e_-e'_)] 



2h J {2ny' 

+cr . (Afe - AfeO [Sie+ - e'_) - S{e. - e'+)]} (13b) 



Jo{Sk) = — ^ j j^^\Ukk'\^{Sk - Sk')5{£(ik - £ok')- (13c) 

The full energies e± = eofc i Ak/2 and e'^ = £ofe' i Afc//2, where = |Afc|. The scalar term Jo (rife) reproduces 
the ordinary Boltzmann-equation scattering term. The side-jump term (rife) constitutes a correction that reflects 
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the presence of the spin-orbit interaction energy with the electric field in the condition for energy conservation. We 
expand the (5-functions in this scattering term in as 

^(e+ - e^,_) = 5{eok - £ofe') + - S{eok - sok'), (14) 

with corresponding expressions for the other combinations of (5-functions. Adding all contributions together the 
scattering term J"- (nfe) simplifies considerably and we obtain the final expression 

■^sl ('^fe) = — ^ |Wfefe'|^(nfe-nfe/)icr-(Afe- Afe/)^^J(£ofe-eofe')- (15) 

I 



The presence of this scattering term reflects the fact that 
the total energy including the spin-orbit interaction en- 
ergy with the electric field is conserved in elastic colli- 
sions. 



B. Scattering correction due to the change in r 
during collisions 

We have so far ignored the presence of the term eE ■ r 
in the time evolution operator. At this stage we would 
like to determine the additional scattering term linear 
in E arising from it, which we denote by J^-^ (ufc). To 
accomplish this we use Eq. (7) to find the correction g^^, 
to gttk' arising from the presence of eE ■ r in the time 
evolution operator. Using Eq. (9a) we will then obtain 



Jl^iuk) as {i/h) [Hu.9%k- 

Starting from Eq. (7), we expand the time evolution 
operator to first order in the term eE ■ r. Using the ma- 
trix elements of the ordinary position operator f between 
Bloch states 

{k\f\k')=i^S{k-k'), (16) 

we obtain additional terms of the form 

t'e--^^^'eE.^Sik-k') = (^e---') eE.-^^5ik-k'). 

(17) 

We integrate over t' as before, and after some lengthy 
but straightforward algebra, we obtain 



We have not written out explicitly a contribution to g^.^, containing terms of the form duk/dk. We find that such 
terms drop out in the final evaluation of spin currents when the scattering potential is elastic, as a result of integrating 
over k and k' . In the final analysis these terms involve the square matrix element \Ukk'\^ which does not have any 
contributions linear in A. We find that the leading contribution due to these terms is oc and may therefore be 
neglected. Evaluating the fe-derivatives of the impurity potentials gives 

Substituting this into the Eq. (18) and subsequently evaluating (nfe) = (i/h) [Hu,g^]kk we obtain the scattering 
term 

J,^ (rife) = ^^^^^ J ■^;^\Uk'k\^{nk-nk')E-crx{k-k')-^—S{eok-sok'), (20) 



which is easily seen to be exactly equal to j^(nfe). The 
sum of these terms constitutes the total side-jump scat- 
tering term Jsj(»^fc) = Jsj{nk)+Jsj{nk) = 2Jg"(nfe), which 
contains the well-known factor of two associated with side 
jump.^^ We emphasize that we obtain this reinforcement 



of the side jump directly from the scattering term, and 
our work shows no evidence that it is related in any direct 
way to the integral of the velocity operator [see Eqs. (29) 
and (35) below] over the time of a coUision.^^ 
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IV. CONTRIBUTION OF THE SIDE-JUMP 
MECHANISM TO THE SPIN-HALL CURRENT 

Wc have derived a contribution linear in the electric 
field to the scattering term appearing in the kinetic equa- 
tion. This contribution is brought about by the spin- 
dependent interaction of the charge carriers with the elec- 
tric field due to the spin-orbit interaction. In this section 
we will first evaluate the correction that this term yields 
in the spin-dependent part of the density matrix, and we 
will show that this correction accounts fully for the side- 
jump spin-Hall current including the important factor of 
two.^^'^°''^'' This is done in subsection IV A. To show that 
this is the only side-jump contribution to the spin-Hall 
efi'ect, subsection IV B will demonstrate that the modi- 
fications to the velocity operator do not contribute any 
additional terms to the spin current. 



A. Contribution of the side-jump scattering term 

We need to find the contribution to the spin-Hall cur- 
rent brought about by the additional scattering term 
Jsj{nk)- In order to further evaluate Eq. (8) we let 
fk = /ofe + fsk, where /ofc(eofc) is the Fermi-Dirac func- 
tion, which is a scalar in spin space in the case un- 
der study, and fsk is the correction we will determine 
from the kinetic equation. Firstly, the kinetic energy 
part of the Hamiltonian, Sofe: drops out of the com- 
mutator. Also, in the commutator [i <t • Afc,/^] we 
note that | <r • Afc is first order in the electric field, so 
the density matrix can be replaced with the equilibrium 
Fermi-Dirac function, which is a scalar in spin space, so 
[| <T- Afc, /ofc] = 0. Moreover, in the side-jump scattering 

term Jsj(?T.fc), which is also first-order in the electric field, 
we may replace rifc by /ofc. Following some short and 
straightforward algebra, the side-jump scattering term 
can be written as 

li{fok) = -T-Ak6is-SF), (21) 

Tp 

where Tp is the usual momentum relaxation time. 

Next, we decompose fsk into a part scalar in spin 
space and a spin-dependent part, thus fsk = J^Bfel + 
Ssk- The equation for nsk in the steady state is the 
ordinary scalar Boltzmann equation. The term [He, f]k 
becomes {e/h)E ■ (dfo/dk) which is a scalar and acts as 
the source term for nsk- We write the scalar part of the 
Boltzmann equation as 



(22) 



where the scattering term Jo{nk) has been defined in Eq. 
(13a). The solution for nsk is written as 



riEk = 



TpcE dfok _ hTp eE ■ k dfpk 
h dk m* deok 



(23) 



The equation for 5*^;;, in the steady state is Jo{SEk) = 
—Jsjifok), in which the RHS, Jsj(/ofc), acts as a source 
term for SEk- Substituting the explicit expressions for 
the two scattering terms, this equation can be written in 
a simpler form as 



= cr • Afc d(£ - sf), 



and has the simple solution 



SEk = -o" • Afc d{e - Sf)- 



(24) 



(25) 



We would like to draw attention to the fact that the ad- 
ditional side-jump collision integral is part of the source 
for Ssk- The source term itself contains a factor of 1/rp, 
which cancels the 1/rp appearing on the LHS of the equa- 
tion for SEk- This explains why the end result for the 
side jump contribution to the spin-Hall current does not 
depend on the strength and shape of the impurity poten- 
tial. 

Now that we have found the solution for SEk, in other 
words the spin-dependent part of the density matrix in 
an electric field, we can determine its contribution to the 
spin-Hall current. We denote the components of the spin 
velocity as which corresponds to a spin component i 
flowing along the direction j. In finding the contribu- 
tion to the spin-Hall current due to SEk we can restrict 
ourselves to the term in the spin velocity to zcroth or- 
der in the electric field, which is Vj = {hkj /m*) hai/2. 
For E = (Ex, 0,0) the side-jump Hamiltonian |cr • Afc 
becomes 



eXE ■ cr X k = —eXExkyUz, 



(26) 



This gives a side-jump spin-Hall current as a result of the 
modification to the scattering term 



(fik Hk 



(27r)" 



^ tr (a, SEk)=neX Ex, (27) 



where the trace is taken over the spin components and 
n is the density. This term therefore gives a spin-Hall 
conductivity <Tyx\^''^ = neX, which is the usual side-jump 
term in the spin-Hall current. ^'''^"■^'''^^'^'''^^ This result is 
valid in both two and three dimensions. 



B. Vanishing contribution of the corrections to the 
velocity operator 

It will be shown in this subsection that the correction 
to the velocity operator linear in the electric field does 
not yield any additional terms in the spin-Hall current. 
The velocity operator is defined as the time derivative of 
the physical position operator, which in turn is given by 



»^phys = r + X& X k. 



(28) 
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Notice that all spin-orbit interactions can be most di- 
rectly derived by replacing r by Tphys in the direct inter- 
actions and expanding to first order in A. 

The velocity operator has a part which is diagonal in 
fe and is given by 

= — + - [eXE ■ & X k,f]k + - [eE ■r,X& x k]k 

m* n h 



hk 

m* 



2eA 



cr X E. 



(29) 



The fc-diagonal part of the spin velocity Vj, up to first 
order in the electric field, is thus 



,• hk-j h<Ji eX , 
v} = ^ — --{i<rxE)„a,} 



(30) 



The disorder potential is the potential U due to the full 
ensemble of impurities present in the system. In the final 
result for the spin-Hall current, the fc-off-diagonal part 
of the velocity operator will be traced with the fc-off- 
diagonal part of the density matrix gkk', which in the 
first Born approximation is also linear in U. Once this 
is done, a configurational average will be performed over 
the impurities. In the end we seek the result to first order 
in A. However, it proves more straightforward to work in 
terms of the full potential U until the end. Only then we 
will restrict ourselves to the terms which are first order 
in A. 

With these insights in mind we proceed, Eq. (33) yields 



Vkk' = -^cr x {k- k')Ukk', (35) 



where {A, B} = AB + BA. For an electric field along x 
we obtain for the fe-diagonal part of the spin velocity Vy 
the expression 



h hky ^ „ 

2 m* 



(31) 



The £^-dependent part of the spin velocity operator is a 
scalar in spin space, and its contribution to the spin-Hall 
current is found by multiplying by the scalar part of the 
equilibrium density matrix /ofc. It gives us the term 



vel.d 



= -eX 



(2^) 



^ tv fok = -neXE^ (32) 



so its contribution to the spin-Hall conductivity is 
a^J-M = _„eA. 

This is, however, not the full story. The velocity op- 
erator also has a term that is off-diagonal in the wave 
vector, which is referred to as Vkk' and is given by 



Vkk' = ^ [Hu , rphyslkk' 



(33) 



The matrix elements of the impurity potential are given 
by Eq. (5). The part of the matrix element Vkk' origi- 
nating from r is easily seen to be 



dUkk' dUkk' 
dk' dk 



(34) 



where we have written the matrix elements of the full 
potential U. Note that this expression has not been av- 
eraged over impurity configurations. This expression for 
the fe-off-diagonal part of the velocity operator holds be- 
cause, for elastic scattering, the scalar part of the scat- 
tering potential Ukk' = W(fc — fc') depends only on the 
difference fc — fc'. Its contribution to the velocity oper- 
ator is immediately seen to be; zciro. This can also be 
understood by noting that the scalar part of the impu- 
rity potential commutes with r and does not contribute 
to the velocity. 

The fc-ofF-diagonal part of the velocity operator con- 
tributes to the side-jump spin-Hall current. To find its 
contribution we return to Eq. (7) and integrate over time 
to find 



gkk' = i'n-S{eok - £ok')Ukk'{fk - fk')- 



(36) 



This expression also contains the matrix elements of the 
full impurity potential and has not been averaged over 
impurity configurations. 

The contribution of the fc-off-diagonal part of the ve- 
locity operator to the spin current is found by taking 
the trace of the spin velocity arising from Eq. (35) with 
the fe-off-diagonal part of the density matrix given in Eq. 
(36). This yields for the spin-Hall current 



vel^od 

= - tia, 
2 



d^k 



(27r) 

d-^k 



d'^k' 



Tvlu,9k'k 



(27r) 



^{kx - k'^) {Ukk'Uk'k) 6{sok' - eofe)(/fc - fk'), 



(37a) 



(37b) 



where the bracket denotes the average over impurity configurations. At this stage we introduce the simplification that 
we require only terms to first order in A. We note that, since the entire term in Eq. (37) already contains A, the other 
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terms in this equation are needed only to zeroth order in A. Consider first the term proportional to fc^, which, to first 
order in A, can be written as 



-iXh) 



27r 



(2^ 



{Ukk'Uk'k) 5{eok' - eok){nk - ny) 



= -{Xh) 

= -m J 



d'^k 



(27r) 
d'^k 



-^kxJo{nk) 



kx 



h dkr 



(38a) 
. (38b) 



where the last replacement follows from the scalar Boltz- 
mann equation, as written in Eq. (22), assuming, as bo- 
fore, that E II X. Further, Eq. (37) also contains a term 
proportional to k'^ , which is easily seen to give exactly the 
same contribution if one swaps k and fe' in the summa- 
tion. The contribution of the fc-off-diagonal part of the 
velocity operator to the spin Hall current to first order 
in A is therefore 



veLod 



= -2XeE^ 



d^k 
(2^ 



Ok 



dkx 



= neXEx. (39) 



The spin-Hall conductivity originating from this term 
neX and it exactly cancels the contribu- 



z \vel^od _ 
yx\ 

z \veLd 



IS a. 

tion CTyj,!"^'''^ from the fc-diagonal £?-dependent part of 
the velocity operator. The physical explanation of this 
cancellation is the following. Notice that the full spin- 
dependent part of the velocity operator, from Eqs. (29) 
and (35) is 



2eA 



cr X E 



2iX 



(T X (fe - k') Ukk' 



(40) 



which contains the total force acting on the system. Ac- 
cording to Ehrenfest's theorem, the expectation values of 
position and momentum obey time evolution equations 
analogous to those of classical mechanics. Consequently 
the expectation value of the force should be zero in the 
steady state, consistent with the earlier suggestion that 
the total force acting on the system does not contribute 
to the spin current.^'' We note also that the presence of 
the velocity terms off-diagonal in wave vector is crucial 
in obtaining the correct side-jump contribution in the di- 
agrammatic Kubo-formula approach, as demonstrated in 
Ref. 20. 



V. SUMMARY AND DISCUSSION 

We have completed the formal derivation of the side- 
jump spin-Hall current, where we have considered (in the 
absence of intrinsic spin precession) all contributions to 
the kinetic equation in the first Born approximation. We 
will now discuss our findings and their implications. It is 
evident from our analysis that, within the kinetic equa- 
tion framework, the side-jump spin-Hall current origi- 
nates solely from the modification of the Boltzmann col- 
lision integral due to the spin-dependent interaction en- 
ergy of an electron with the external electric field and 



the impurity field. The density-matrix formulation of 
the problem shows that the spin-orbit interaction with 
the electric field alters the condition for energy conser- 
vation, since the total energy conserved during collisions 
must include the spin-dependent part. In addition, the 
spin-orbit coupling with the impurities causes a change 
in the position of the electron during scattering process, 
which again affects the energy balance via the interac- 
tion energy eE ■ r. This effect doubles the size of the 
side-jump current. 

The understanding of the side-jump effect emerging 
from this derivation differs from the conventional expla- 
nation, according to which this phenomenon is attributed 
to the linear-in-£J modification of the velocity operator. 
It is clearly seen in the previous section that the full ve- 
locity operator in the presence of spin-orbit interactions 
contains an extra term due to the impurity potential, 
which is off-diagonal in wave vector. This term, and the 
fc-diagonal velocity operator result in two corrections to 
the spin current that are equal in magnitude but op- 
posite sign so they cancel out. This could be justified 
informally using the fact that the velocity operator con- 
tains the total force acting on the system (or rather the 
term to leading order in A of this force), and therefore 
should vanish in the steady state^'^. Therefore, contrary 
to conventional assumptions, the side-jump contribution 
to the spin-Hall current is traced to the qualitatively dif- 
ferent carrier spin dynamics during collisions. It is not 
traced to any Xme&r-'m-E correction to the velocity oper- 
ator. The linear-in-£J correction to the velocity operator 
is canceled by the contribution of the off-diagonal in fe 
correction to the velocity operator. The importance of 
this off-diagonal-in-fc correction in the velocity operator 
in obtaining the correct side-jump current was recognized 
in the diagrammatic Kubo-formula approach as shown in 
Ref. 20. 

We emphasize that in obtaining our results we have 
used a rigorous quantum mechanical formulation, start- 
ing with the quantum Liouville equation and making, 
in the course of the derivation, the same assumptions 
that are characteristically made in linear response theory. 
Therefore our formalism could in some sense be regarded 
as being built from the ground up. The approach we have 
used demonstrates that the derivation of the side jump 
does not need to rely on intuitive semiclassical ideas as 
long as the collision integral is derived rigorously from 
the fundamental starting point of all transport theories. 
In contrast, Ref. 37 using semiclassical Boltzmann argu- 
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mcnts counted six terms contributing to side jump, but 
did not clearly indicate which terms should cancel. Since 
for the conduction band the side-jump contributions have 
the same magnitude but opposite signs, it led to freedom 
in choosing which terms cancel and still obtaining the 
correct amplitude of the side-jump. For example Ref. 19 
counted the terms from the anomalous velocity and the 
shift of the position operator. Our analysis shows exactly 
which terms arc non-zero and which contributions can- 
cel and is in agreement with the Kubo derivation of side 
jump contributions presented in Refs. 20,49. Our results 
are also in agreement with recent extensive studies of the 
anomalous Hall effect in magnetic semiconductors. ^^'^^'^-"^ 
In particular, Refs. 42 and 52, by comparing the semi- 



classical description of the side-jump with the Kubo and 
Keldysh formalisms, found that the velocity operator is 
unchanged and that the two contributions to the side 
jump have different origins due to the renormalization of 
the distribution function and due to the modification of 
the conservation of energy. 
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